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Abstract 

The quantum mechanical formalism for position and momentum of a par- 
ticle in a one dimensional cyclic lattice is constructively developed. Some 
mathematical features characteristic of the finite dimensional Hilbert space 
are compared with the infinite dimensional case. The construction of an un- 
biased basis for state determination is discussed. 

I. INTRODUCTION 

In the quantum mechanical description of physical systems, it is often assumed a contin- 
uous set of states requiring an infinite dimensional Hilbert space. There are however many 
physical systems whose states belong to a discrete and finite set, with a quantum mechanical 
description formalized in a finite dimensional Hilbert space. The best known example of this 
is the quantum treatment of angular momentum, indeed, a paradigm for quantum mechan- 
ics, presented in every text book. A less known example of discrete quantum mechanics, 
presented in this work, involves the description of position and momentum observables in a 
finite dimensional Hilbert space. 

In this case, the position observable does not take values in a continuous set but instead 
it can take values on a lattice. One important practical motivation for studying such systems 
is that any computer simulation of position and momentum will necessarily involve a finite 
number of sites. On the other side, a highly speculative motivation is that the posible 
existence of a fundamental length scale, that is, a measure of length below which the concepts 
of distance and localization become meaningless, can make a discrete quantum mechanics 
more appropriated than a continuous one. 
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In order to make this work useful from the didactic point of view, the formalism of 
quantum mechanics in a finite dimensional Hilbert space will be presented in a constructive 
way where all steps are logically connected. This work may therefore be a useful comple- 
ment to any text book where quantum mechanics in infinite dimensional Hilbert space is 
developed. Finally, another didactic feature of this work is that finite dimensional quantum 
mechanics requires many interesting mathematical tools such as some finite sums and the 
Discrete Fourier Transform that are not usually presented at the undergraduate level. Fur- 
thermore, the important differences between finite and infinite dimensional Hilbert spaces 
are emphasized and the limit when the dimension becomes infinite is considered. 

II. NOTATION AND DEFINITIONS 

In this work we will consider a particle in a one dimensional periodic lattice with N sites 
and lattice constant a. The quantum mechanical treatment of this system requires an N 
dimensional Hilbert space 7i. Given any two elements of this space $ and \1/ we will denote 
their internal product by ($, ^). We will use operators of the form A = •), where the 
dot indicates a space holder to be occupied by the Hilbert space element upon which the 
operator acts. The corresponding hermitian conjugate is A' = ■). 

Although we don't need to choose any particular representation for the abstract Hilbert 
space H,, it may be convenient, for didactic reasons, to specialize the formalism in a three or 
four dimensional Hilbert space whose elements \1/ are column vectors of complex numbers. In 
this case •) represents a complex conjugate row vector and operators are square matrices. 
This special representation is recommended for clarity but it is important to emphasize to 
the students that the formalism of quantum mechanics can be construed in the abstract 
Hilbert space and a particular representation is never required. The mathematical beauty 
of quantum mechanics is, indeed, most apparent in the abstract formulation. 

Any basis {<Pk} in the Hilbert space will have N elements labeled by an index k running 
through the values —j, —j + 1, — j + 2, • • • , j — 1, j, with j = |, 1, |, 2, |, • • • corresponding 
to the dimensions N = 2j + 1 — 2, 3, 4, 5 • • •. This choice of labels has some advantage and 
some disadvantage. The main virtue of this symmetric labeling is that it corresponds to 
the physical concepts of position and momentum that can take positive and negative values. 
The main shortcoming of it, is that there are many summations and results that are usually 
given in books with integer indices running from to N — 1. In order to cure this deficiency, 
we give in Appendix A some of these sums with the symmetric index. Furthermore, with 
this notation we must keep in mind that for even N, j takes half-odd-integer values and this 
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may be relevant in modular mathematics. 

We will adopt a very useful notation for the principal N th root of the identity defined by 



; 2TT 



lu = e 1 ^ . (1) 
Integer powers of this quantity build a cyclic group with the important property 

l = UJ Nn = ujM +1 > , Vn = 0, ±1, ±2, • • • . (2) 



III. POSITION AND MOMENTUM 

The position of the particle in the lattice can take any value (eigenvalue) ax where a has 
units of length and the discrete number x can take any value in the set {—j, — j + 1, • • • , j — 
The state of the particle in each position is represented by a Hilbert space element ip x 
and the set {<p x } is a basis in 7i. In the spectral decomposition, we can write the position 
operator X as 

j 

X = J2 aX( Px{<Px: •) , (3) 
x=-j 

that clearly satisfies Xip x = axip x . We can now construct a translation operator T with the 
property 

T = I Vx+1 ,X ^ J (A) 

We will later explain the reason for defining this operator periodic for odd dimension but 
antiperiodic for even dimension. This operator is given by 

T= Vx+iiVxr) + (-ly^V-i^jv) , (5) 

x=-j 

with its hermitian conjugate 

Tt = £ <P*(<Px+i, •) + (-if-^AV-j, •) • (6) 

x=-j 

It is straightforward to check that this operator is unitary, TT^ = T^T = 1, and therefore its 
eigenvalues are complex numbers of unit modulus and their eigenvectors build a basis (see 
Appendix B). Let then {4> p } and {A p } for p = — j, — j + 1, • • • , j — 1, j, be the eigenvectors and 
eigenvalues of T. In order to determine them, we expand P in terms of {<p x } and consider 
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t E ifPxApiVx = \Y1 {(pxAp)?* ■ ( 7 ) 

x=-j x=-j 

From Eq.([|) we get 

{<Px-i, <t> P ) = *p(<Px, <P P ) for x ^ j , ^ 

(<Pj,<f> P ) =\{ ( P-jA P ) ■ 

Up to an arbitrary phase, that can be absorbed in the definition of <f) p , and considering the 
normalization of (f) p , the solution of the above equations is 

((fx, <P P ) = -7= uj px , and X p = u~ p . (9) 
v Jy 



The two bases {(f x } and {4> p } are then related by 



ft = it • ( 10 ) 



1 E . (ii) 



and 



N ~ 

Except for the symmetric index and a different factor, this is essentially the Discrete Fourier 
Transform. Notice that if we had not defined the translation operator antisymmetric when 
iV is even, then we would not have obtained such a simple relation in Eq.@ above and we 
would have obtained different expressions for N even or odd. In other words, we choose to 
define the translation operator in a way to obtain a simple relation between the bases. The 
complication in the definition of the translation operator is also related to our choice of using 
symmetric indices. If we had chosen indices running from to N — 1, then a translation 
operator periodic for all N would have lead to the two bases related by the Discrete Fourier 
Transform (also expressed in terms of asymmetric indices). 
Therefore we have 

T<j) p = uj- p <j) p , (12) 

or equivalent, 

T= J2 u~ p M<Pp,') ■ (13) 
p=-j 

We can now construct an hermitian operator P as a superposition of projectors in the basis 
{<f> P } 
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p= E grtpitp,-) , (14) 
p=-j 

where g is a real constant to be determined later. Clearly, this operator is hermitian and 
satisfies the eigenvalue equation P(p p = gp<p p . From this equation and from Eq.(fL2|), and 
doing the power expansion of the exponential we prove that 

T = exp (-43 . (15) 

We can now assign a physical interpretation to the operator P. This last Eq. (|T5|) , together 



with Eq.(^j) means that P is the generator of translations in the position observable. We 
identify this observable P, as is done in classical mechanics, with the momentum. If the 
position observable takes value in a lattice with lattice constant a, then, the momentum 
observable must also assume values in a lattice with lattice constant g. In the next section 
we will see that these values must be related by ga = 2it/N, that is, the momentum lattice 
is the reciprocal lattice of the position. 

In an identical manner as was done before, we can now construct a unitary operator B 
that "boosts" the momentum states 

B4> P ={ ' P ^' , (16) 



and show that 



and 



and also that 



B<p x = cu x ip x , (17) 



B= ^V*<^,-> , (18) 

x=-j 



B = «p(4f). (19) 



IV. COMMUTATION RELATION AND THE LIMIT N -»• oo 

Every quantum mechanics textbook emphasizes that the position and the momentum 
operators satisfy the commutation relation [X, P] = i (we use units such that h — 1). 
However, in most cases it is not mentioned that this commutation relation is false in a 
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finite dimensional Hilbert space. It becomes clear that this must be so because one can 
prove that the commutation relation [X, P] = i implies that the operators X and P are 
unbound. However, in a finite dimensional Hilbert space all operators are bounded; therefore 
such a commutation relation is imposible in a finite dimensional Hilbert space. An explicit 
calculation of the commutator in the position representation, that is, in terms of the basis 
{tp x } results in 

3 3 3 1 / Oir \ 

[X,P]=agY: E E k(s-r)—exp (i—k(s-r))ip s (ip r ,-) . (20) 

k=-j s=-j r=-j iV \ iV / 

The sum over k can be performed but it is advantageous to leave it unperformed. We can now 
see that in the continuous limit, where N — * oo, a — > and g — > but agN — ^constant, the 
above commutator approaches the value i, provided that lattice constants satisfy agN = 2tt. 
In this limit, the sums over discrete indices k, s, r, become integrals over continuous variables 
k, a, p, according to the scheme 

2tt, /2tt /2tt 

N k ^*> v^ s ^ a ' ix r ^ p 

I2n /27 , s 4^ 

f« ^ ' v ^ ~* vip) '5 ^ J-oo ■ 

The continuous limit is then given by 

agN f°° f°° —1 f°° 
[X,P] — da dp— dK iK(a - p)exp(iK(a - p))ip(a)(ip(p),-} . (21) 

Z7T J-oo J-oo Z7T J-oo 

The sum over k, that was left unperformed, assumes in the continuous limit a simple form. 
Indeed, the integral over k is a well known representation of Dirac's delta function. Therefore 

.agN 




[X,P]^i^- Hda fdp S(a - pMa)(<p(p), 

Z7T J-oo J-oo 



= ^/>W( f W..) = ifl, (22) 

where we have used the completeness relation. Therefore the usual commutation relation 
for the continuous case is recovered, provided that 

agN = 2tt . (23) 



V. STATE AND TIME EVOLUTION 

At any instant of time, the state of the particle will be determined by a Hilbert space 
element We can represent this state in the position or momentum representation, that 
is, expanded in the bases {<p x } or {4> p }. 
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* = J2 c *Vx = E d v$v » ( 24 ) 



P=-J 



where the complex coefficients c x and d p are related by the Discrete Fourier Transformation 
rf P = E w ^ > c - = -ffi E > ( 25 ) 



and their absolute values squared \c x \ 2 and \d p \ 2 represent the probability distributions for 
position and for momentum. Let ^(t ) be the state of the system at some instant to, that 
we can choose to be to = 0. In Schrodinger's picture, this state will evolve according to the 
time evolution unitary operator given in terms of the hamiltonian H as 

U t = exp(-zfft) . (26) 

This description of the time evolution is equivalent to Schrodinger's equation if time is repre- 
sented by a continuos variable. However in some cases it may be convenient to assume that 
also time takes discrete values giving preference to the formulation with the time evolution 
operator above. If the state is given in the position or in the momentum representation, 
the coefficients of Eq.([24]) will become functions of time c x (t) and d p (t). Let us consider 
for instance the case of a free particle with hamiltonian H = P 2 /2m. In the momentum 
representation the coefficients are simply given by 

d p (t) = dp(0) exp (^-i^tj = d p (0)uj~ p ^ , (27) 

where we have introduced a time scale t defined by 

- — • (28) 

whereas in the position representation we have 

c r (t) = t *(0)-^ E > {M) - p2 ^ • (29) 

x=-j P=-j 

The second summation in Eq.(^) is a Discrete Fourier Transform that becomes, in the 
continuous limit, a Fourier Integral Transform of a gaussian function with a very well known 
result. In our discrete case, this summation can not be evaluated in general. This is an 
example of the difficulties encountered in the discrete case. It took Gaufi four years working 
"with all efforts"§ in order to evaluate a similar summation (the so called "Gaufi sum" ) for 
some special values of the parameters involved. In any case we can see that the state is 
periodic, + T) = ^(t) with period T = Nt if N es odd and T = 4Nt if iV es even. 
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VI. STATE DETERMINATION AND UNBIASED BASES 



At an early stage in the development of quantum mechanics, Paul! raised the question 
whether the knowledge of the probability density functions for position and momentum 
where sufficient in order to determine the state of a particle. Since position and momentum 
are all the (classically) independent variables of the system, it was, erroneously, guessed 
that this Pauli problem could have an affirmative answer. Indeed, many examples of Pauli 
partners, that is, different states with identical probability distributions, where found. A 
review of theses issues, with references to the original papers can be found in refsJH^. 
Considering the similar problem in classical statistic, we should not be surprised to find 
out that the Pauli question can not have a positive answer. The marginal probability 
distribution functions of two random variables uniquely determine the combined distribution 
function only in the case when they are uncorrelated, that is, when they are independent 
random variables. Position and momentum are, however, always correlated; that is indeed 
the essence of Heisenberg's uncertainty principle, and therefore we should not expect that 
their distributions uniquely determine the quantum state. 

Explicitly stated in our case, the Pauli question is: can we find the set of N complex 
numbers {c x } that determine the state in Eq.(E3) with the knowledge of the sets { | | 2 } 



and {| dp| } related by Eq.(|25f)? Let us notice that the state has an arbitrary phase and is 
normalized; therefore we only need to find 2N — 2 real numbers in order to determine the 
state. The known numbers {|ca;| 2 } and {|d p | 2 } are not independent because the numbers of 
each set are probabilities and they should add to 1. We have therefore 2N — 2 equations 
at our disposal in order to find 2N — 2 unknown. However, the equations are not linear 
and they are not sufficient for an unambiguous determination of the state. There is another 
very important feature in these equations. We will see that not every set of data {|c x | 2 } and 
{|dp| 2 } are compatible. The equations have solution only if the position and momentum data 
satisfy a number of relations. These constraint on the data is just Hiesenberg's uncertainty 
principle and is a consequence of the relations in Eq . (|25|) . These concepts are clarified by 
an example with N = 2. 

Let {ip-,ip + } and {0_,0 + } be the position and momentum bases in two dimensional 
Hilbert space. Their internal product is given by Eq.([|). An arbitrary state, normalized 
and with a phase fixed is determined by 2N — 2 = 2 numbers < g < 1 and < a < 2ir: 



i> = ge ia <p„ + £?V+ • (30) 
The independent data on position is | | 2 = g 2 , that directly determines g and the 
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independent data on momentum is \{4>-,ip)\ 2 = w 2 . With this last data we must determine 
a. Using that (0_,</?±) = exp(±z7r/4)/v2, we get after some algebra that 

^ 2 -i/2 

sin a = — . . (31) 



This equation can only have solution if 

zu 2 - 1/2 



< 1 , (32) 



that, after squaring and arranging, becomes 

(w 2 -l/2) 2 + (g 2 -l/2) 2 < (1/2) 2 . (33) 

This relation is indeed the uncertainty principle: if g 2 = or 1, that is, exact localization 
in ip + or <£>_, then w 2 = 1/2, that is, maximal spread in momentum and, vise versa, exact 
momentum {vj 2 = or 1) implies maximal spread in position (g 2 = 1/2). Now, even if 
the data is consistent with the uncertainty principle, there is an ambiguity in the solution 
of Eq.([H]) because if a is a solution then 7r — a is also a solution. This ambiguity can not 
be solved with the given data and requires more experimental information. We will next 
consider what observables can we measure in order to determine the state without ambiguity. 

From previous example it is clear that we need further information besides the distribu- 
tion of position and of momentum in order to determine the state of the particle. This will 
involve an observable depending on both, position and momentum because any observable 
depending on only one of them will not bring new independent information. Some can- 
didates may be X + P or the correlation XP + PX or any function F(X, P) symmetric 
or antisymmetric under the exchange X <-> P. Perhaps the best choice of an observable 
that provides information on the system not available in the knowledge of position and mo- 
mentum distributions is an observable whose associated basis is unbiased to the position 
and to the momentum bases. Two bases in a Hilbert space are unbiased when they are as 
different as posible in the sense that any element of one basis has the same "projection" on 
all elements of the other basis. More precisely, the modulus of their internal product is a 
constant for all pairs. Unbiased bases are candidates for the quantum mechanical descrip- 
tion of classically independent variables like position and momentum; indeed we have from 
Eq. (|9|) | ((p x , (p p ) | = 1/ vN Vx, p. This leads us to the search of a basis {r] s } unbiased to {(f x } 
and to {4> p }. 

The importance of unbiased bases associated to non commuting observables was recog- 
nized by SchwingeiH long ago but the existence and explicit construction of maximal sets of 
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mutually unbiased bases for any dimension is still an open problem. When the dimension 
iV is a prime number, a set of iV + 1 mutually unbiased bases was presented^! and this was 
extended to the case when N is a power of a prime number!. In order to find unbiased bases 
we will follow the method given by Bandyiopadhyay et al.0. We have seen that the eigen- 
vectors {4> p } of the operator T that produces a translation or shift on the basis {<p x } build 
an unbiased basis to {tp x }- This result is generalized in reference |10| where it is shown that, 
if N is prime, the eigenvectors of the unitary operators T, B, TB, TB 2 , ■ ■ ■ TB N ~ l build a set 
of N + 1 mutually unbiased bases where T and B are the translation operators for position 
and momentum defined in Eqs.@ and (|16|). In our case we want to find a set of only three 
unbiased bases and therefore we just consider the first three operators that provide unbiased 
bases for any iV (prime or not). The first two operators provide the bases and {(p x }, 
that are related by the Discrete Fourier Transform and are clearly unbiased. One can easily 
prove, with Eqs.([|, |17|) and (|16| , |i2|) that the operator TB is a shift operator for both bases 



{4> p } and {ip x } and therefore its eigenvectors {r] s } build a basis unbiased to both of them. 
We have indeed 

TB<p x ={ _ l , (34) 



and 



TB<^=\ ; P+1 ' PTJ , . (35) 

The eigenvectors of TB in the position representation are found by expanding r] s in the basis 
{ip x } and using Eqs.([£|) and the relation 

TB Vs = uj% . (36) 

In this calculation one must use with care the modular mathematics defined in Eq.(§). This 
results in 



1 



j 



Vs = -4= E c> 2 " (s+ ^ . (37) 



N x= ^ 



With a similar calculation we obtain the eigenvectors of TB in momentum representation 

ris = -L E c-^-(s+§)^ . (38) 
VN p= _j 

Clearly we see that {?7 S } is unbiased to {ip x } and to {4> p }- The analytical calculation of 
Discrete Fourier Transforms is much more difficult than the Fourier Integral Transform and 
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therefore it is of mathematical interest that, from the last two equations we can obtain the 
Discrete Fourier Transform for a family of sequences. If we equate Eqs. (p7[) and (^) and 
we expand (f x in terms of {<ft p } we obtain 

1 >A i x 2_ bx _ vx _ i„2_ bp ( b = 0, ±1, ±2, ■ • • , N even 
— = > u;2 u px = u 2 p op ,< . (39) 

\fN x t±j { b = ±1/2, ±3/2, - ■ • ,iVodd 

We rewrite this result in terms of the asymmetric indices more usual in the mathematical 
literature as, 

i N-l 

(-l) n uj^ n2 ~ bn uj- mn = (-l) m uj-^ m2 - bm , V6 integer . (40) 



N n=0 

The unbiased basis found is of course not unique. It was generated by the operator TB 
but there are many other operators whose eigenvectors build unbiased basis to {ip x } and to 
{4> p }- Indeed, any combination T n B m or B m T n where n and m are not divisors of iV could 
do the job. 

We will now try to find the physical meaning of the basis {i] s }. That is, we would 
like to find an hermitian operator S(X, P) with {r] s } as eigenvectors corresponding to the 
eigenvalues hs with h = ag. That is, 

TB = exp(iS(X,P)) . (41) 

In terms of the operators X and P, and using the relation agN = 2n, we have 

exp {iS{X, P)) = exp (-iaP) exp (igX) . (42) 

Notice that here we can not use the Baker-Campbell-Hausdorff relation e p e x = e p+x ~ 1 ^ 2 
that is valid in the N — > oo case, where the commutator [X, P] is a constant. If this 
were possible, then the operator S would be simply equal to gX — aP. This is one of the 
subtle differences of finite and infinite dimensional Hilbert spaces. It is posible to find the 
eigenvectors of the operator X — P but the basis so obtained is not unbiased to either {(p x } 
nor {4> p } however it becomes unbiased in the infinite dimensional limit0. The relation of 
S to X and P is not simple but we can prove that S is antisymmetric under the exchange 



P <->• X and a <-> g. Indeed, from the hermitian conjugate of Eq.(f4~2"D we have 



exp (—iS(gX, aP)) = exp (—igX) exp {iaP) = exp (iS(aP, gX)) . (43) 
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VII. CONCLUSION 



In this work we have presented the quantum mechanical formalism for position and 
momentum of a particle in a one dimensional cyclic lattice in a way that may be useful 
for a didactic complement of the infinite dimensional case presented in quantum mechanics 
text books. In doing this, several mathematical subtleties related to the difference between 
infinite and finite dimensional Hilbert spaces, and of modular mathematics, arise. We have 
discussed the physical and mathematical relevance of unbiased bases and, as consequence 
from the construction of such a basis, the Discrete Fourier Transform for a family of sequences 
is given. 

It is a strongly recommended exercise to reproduce all this work in terms of the asymmet- 
ric indices running from to N — 1. One can see thereby the need to define the translation 
operator always cyclic in order to get the position and momentum bases related by the Dis- 
crete Fourier Transform. The calculations of the eigenvectors of the operator TB are useful 
exercises for the modular mathematics. 

This work received partial support from "Consejo Nacional de Investigaciones Cientificas 
y Tecnicas" (CONICET), Argentina. 



VIII. APPENDIX A 



All sums appearing in this appendix can be derived from the fundamental expression 



N-l 

E- fe 

fc=0 



N 



l-Z 



(44) 



for any complex number z, that in terms of the symmetric index becomes 



k=-j 



z N/2 _ z -N/2 

z 1 ' 2 - z- 1 / 2 



for < 



1=1 l3o 

J 2' ' 2' ' 



N = 2j + l = 2,3,4,5, 
z complex 



(45) 



If z takes the special values z = exp(i-^r) = lo t with r an arbitrary number, then 



E 

k=-j 



^kr = sin(Trr) 



sin(f r) ' 



(46) 



In our case, the number r will often assume integer or half-odd-integer values. For these 
cases we have, 
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_l)»(Ar-i)jV for r = n N ) n = 0, ±1, ±2, • • • 
^cj fcr =<jo for r = ±1, ±2, ±3, • • • ^ nN . (47) 



forr = ±±±i,±-? 



1 _ Ljr ^ _^ 2 , _^ 2 , _^ 2 , 

The first two cases correspond, with the asymmetric index, to the well known result 
J2k=o ukr = N5 r:n N where we see that this choice leads to simpler mathematics. The third 
case above must be handled with care in numerical evaluations because the numerator is 
the fourth root of exp(i^u), with u an odd integer. Therefore it has four posible numerical 
results. The denominator has also two posible results. A formal derivative of Eq.(^6|) with 
respect to the parameter r leads to the result 

V- , kr i sin(Trr) cos(f r) - Acos(Trr) sin(f r) 

Deriving again with respect to r, we can obtain other summations involving higher powers 
of k. 



IX. APPENDIX B 



In most textbooks it is proven that the non degenerate eigenvalues of an hermitian 
operator are real and their eigenvectors are orthogonal. We give here the corresponding 
proof for unitary operators. 

Let T be an unitary operator and A& and 0^ the non degenerate eigenvalues and normal- 
ized eigenvectors. Then, we will prove that, |Afc| 2 = 1 and (fa, fa) = 5 r k- 
From Tfa = Xk4>k and T^T = 1 it follows that 

|A fc | 2 = (Tfa, Tfa) = (fa, T^Tfa) = 1 . (49) 

In order to prove the orthogonality consider that 

Tfa = X k fa -> (fa, Tfa) = X k (fa,fa) , 
T ] fa = X* r fa —> (T ] fa,fa) = X r (fa,fa) . 

Subtracting both equations we get = (A& — X r )(fa,fa). Since the eigenvalues are non 
degenerate, the product (fa., fa) must vanish for k ^ r. 

Since T is bounded, the completeness of the eigenvectors can be proved in the usual way 
and therefore {4>k\ is a basis. 
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